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Abstract 

We compute the decoupling constant Cm relating light quark masses of effective 
n;-flavour QCD to (n^ + l)-flavour QCD to four-loop order. Immediate applications 
are the evaluation of the MS charm quark mass with five active ffavours and the 
bottom quark mass at the scale of the top quark or even at GUT scales. With the 
help of a low-energy theorem Cm can be used to obtain the effective coupling of a 
Higgs boson to light quarks with five-loop accuracy. We brieffy discuss the influence 
on T{H ^ hi). 

PACS numbers: 12.38.-t, 12.38.Bx, 14.65.Dw, 14.65.Fy 


1 Introduction and notation 

Perturbative calculations in QCD are quite advanced and have reached, at least for some 
observables, the four and even five-loop level (see Refs. DEI for a recent review). This 
concerns in particular the renormalization group functions which have been computed at 
four loops in Refs. [3HZ1. The first five-loop result has been obtained recently in Ref. [8] 
where the quark mass anomalous dimension has been computed to this order. 

In order to consistently relate the quark masses and strong coupling constant evaluated at 
different energy scales, both the renormalization group functions and also the decoupling 
relations have to be available. The latter take care of integrating out heavy quark fields. In 
fact, V-loop running goes along with (iV — l)-loop decoupling. Thus, besides the five-loop 
anomalous dimensions also the four-loop decoupling relations are needed. In Refs. PITO] a 
first step has been undertaken in this direction and the four-loop decoupling constant for 
as has been computed (although the five-loop beta function is not yet available). In this 
paper we complement the result by computing the four-loop corrections to the decoupling 
constant for the light quark masses, which supplements the five-loop result for 7m [8]. 


In Ref. m a formalism has been derived which allows for an effective calculation of the N- 
loop decoupling constants with the help of A^-loop vacuum integrals. In the following we 
present the formulae which are relevant for the calculation of the quark mass decoupling 
constant. 

The bare decoupling constant is defined via the relation 




( 1 ) 


where and m®' are the bare quark mass parameters in the full nf- and effective ni(= 
Uf — l)-flavour theory. Introducing the renormalization constants in both theories leads 
to the equation 





= Cmmgifi), 


( 2 ) 


which relates finite quantities and defines (m- Note that primed quantities depend on 
and non-primed quantities on ■ Four-loop results for Zm and Z'^ can be found 
in Refs. ME] and Cm can be computed with the help of 


1 - 

1 + s"^(0) ’ 


(3) 


where S^^(O) and Sy^(O) are the scalar and vector parts of the light-quark self energy 
evaluated at zero external momentum. The superscript “h” reminds that one has to 
consider only the hard part which involves at least one propagator of the heavy quark. 

In the next Section we discuss the calculation of Cm ^^cid its renormalization to arrive at 
Cm- Section [3] applies a low-energy theorem to derive, from the four-loop result of Cm, the 
effective Higgs-fermion coupling constant to five-loop order. We summarize our hndings 
in Section m 


2 Decoupling for light quark masses 

In this section, we compute the decoupling constant Cm ^^id combine it with the four- 
loop result for Zm to obtain the finite quantity Cm- The computation of Cm requires 
the knowledge of the hard contribution to the scalar and vector part of the light-quark 
propagator, see Fig. [Tjfor sample Feynman diagrams. The first non-vanishing contribution 
arises at two loops where one diagram contributes. At three-loop order there are 25 and 
at four loops we have 765 Feynman diagrams. 

The perturbative expansion of Eq. ([3D to four loops leads to 

c = + + , (4) 

where in the last term on the right-hand side only two-loop expressions for and 

S^(0) have to be inserted. 
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Figure 1: Sample Feynman diagrams contributing to the hard part of the light-quark 
propagator up to four loops. Solid and curly lines denote quarks and gluons, respectively. 
At least one of the closed fermion loops needs to be the heavy quark. 


We generate the Feynman diagrams with the help of QGRAF [12]. FORM [T^IT^ code is then 
generated by passing the output via q2e [I^ITG] . which transforms Feynman diagrams into 
Feynman amplitudes, to exp [T5l[T6]. After processing the latter one obtains the result as 
a linear combination of scalar functions which have a one-to-one relation to the underlying 
topology of the diagram. The functions contain the exponents of the involved propagators 
as arguments. At this point one has a large number of different functions. Thus, in the 
next step one passes them to a program which implements the Laporta algorithm HZ] 
and performs a reduction to a small number of so-called master integrals. We use, for 
the latter step, the C++ program FIRE [TB]. Our four-loop result is expressed in terms of 
13 master integrals which we take from Ref. [IB] (see also [2U1 - I22] and references therein). 
All e coefficients are known analytically in the literature except the term of integral 
Jq ^2 (in the notation from Ref. [19]) which has been provided from [23] . 

Note that for our calculation we have used a general gauge parameter ^ of the gluon 
propagator. At four loops, in intermediate steps terms up to order are present, however, 
in the hnal result for all ^ terms drop out. The last term on the right-hand side of 
Eq. dlj) is separately ^-independent since at two loops S^'^(O) and Sy^(O) are individually 
^-independent. The results up to three-loop order have been checked with the help of 
MATAD [21] which avoids the use of the program FIRE since it implements the explicit 
solution of the recurrence relations. 

To obtain we have to renormalize and the heavy quark mass rrih to two-loop order. 
The corresponding MS counterterms are well-known (see, e.g.. Ref. 0). Cm still contains 
poles in e which are removed by multiplying with the factor Z^/Z'^ (see, Eq. ([2|)) which is 
needed to four-loop order mw- Note that Z'^ depends on the strong coupling constant 
of the effective theory, whereas Z^ and Cm expressed in terms of In 

order to achieve the cancellation of the e poles the same coupling constant has to be used 
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in all three quantities. We have decided to replace in favour of which is done 

using the corresponding decoupling constant (as three-loop order mi- Note, however, 
that higher order terms in e are also needed since (as gets multiplied by poles present 
in Z'^. Up to two-loop order they can be found in Refs. |25|,[26]; the three-loop terms of 
order e can be extracted from Refs. 

Our hnal result for the decoupling constant parametrized in terms of the MS heavy quark 
mass, iJih = rrihifi), reads 
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a^r^\mh) 
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(5) 


TT 


(6.850 - lA66ni + 0.05616nf) , 


with In the analytic expression (^(n) denotes the Riemann zeta fnnction 

evalnated at n and a„ = Li„(l/2). 

Often it is convenient to express Cm in terms of the on-shell heavy qnark mass, M^- The 
corresponding analytic expressions are obtained from Eq. ([5]) with the help of the two-loop 
relation between mh{fi) and Mh which can be fonnd in Refs. [27l - l2^ . We refrain from 
showing the corresponding analytic result and restrict the presentation to the numerical 
expression which is given by 


yos = 1 + 

Sm ' 


+ 


a. 


O/)' 


a. 


71 

O/)' 


TT 


2 2 

0.2060 - 0.1389 In iU + 0.08333 In^ A; 

Ml Ml 


2 2 2 

1.477 - 0.9550 In ^ + 0.7384 In^^ + 0.1343 ln=^^ 
Mi Mi Mi 


ni 


2 2 
0.02473 - 0.1227 In - 0.009259 ln=^ 

Ml Ml 


+ 


a. 


Af)' 


TT 


0.2233 


2 2 2 2 
2.674 In ^ + 6.227 In^ ^ + 2.165 ln=^ ^ + 0.2648 In^ ^ 
Ml Ml Ml Ml 


ni 


2 2 
-1.504 -0.6470 In0.9260 In^ ^ 
Mi 


Ml 


- 0.03472 In^ 


0.001157 Ih 


Ml 

4 

Ml 


+ nf 


- 0.1632 In' 


/i 


3 + 

Ml 


0.05616 + 0.005016 In ^ + 0.02392 In 

Mi 


2 + 

Ml 


( 6 ) 


On the webpage [30] we provide analytic results in computer-readable form for a general 
SU{Nc) gauge group. 

In the remaining part of this section we discuss two applications which involve the evalua¬ 
tion of light quark masses at high scales. In the hrst one we compute the running bottom 
quark mass at the scale + = Mt, where Mt is the top quark pole mass. mb{Mt) appears as 
an intermediate step in analyses concerned with Yukawa coupling unihcation. Here the 
role of the heavy quark is taken over by the top quark. In the second application we cross 
the bottom threshold and evaluate the charm quark mass for /i = Mz using mi'^^(3 GeV) 
as input. As input parameters for the numerical analyses we use 

af\Mz) = 0.1185, 


mf^ ) = 
mW(3GeV) = 


5 


0.1185, 
4.163 GeV, 
0.986 GeV. 


(7) 










As a first phenomenological application we consider the evaluation of the bottom quark 
mass at the scale of the top quark with six active flavours using rnf\rnf'^) as input. We 
are interested in the dependence of on the decoupling scale of the top quark. 

Since this scale is unphysical it should get weaker after including higher order corrections. 
Our results, which are shown in Fig. [2al are obtained using the following scheme, where 
N G {1, 2, 3,4, 5} refers to the number of loops: 

• Use Wloop running: —?■ rrif\jjL^^‘^) 

• Use {N — l)-loop decoupling: 

• Use Wloop running mf\Mt) 

The values for as involved in this procedure, and 

are obtained from ai^\Mz) using the same loop-order for the running and 
decoupling as described above for the bottom quark mass. 

In Fig. |2a]m[®^ (Mt) is shown as a function of the scale where the transition from five- to 
six-flavour QCD is performed normalized to the on-shell top quark mass. For the on-shell 
top quark mass we choose Mt = 173.34 GeV [33]. We vary by a factor of 10 around 

the central scale = 1. The one-loop result leads to ^ 2.9 GeV and is not 

shown in the plot. One observes that already the result where two-loop running is used 
(short-dashed line) shows only a weak dependence on pf®®. It becomes even weaker at 
three and four loops (results with higher perturbative order have longer dashes) and results 
in an almost flat curve at five loops (solid line) which can barely be distinguished from 
the four-loop curve. The five-loop results depends on the unknown five-loop coefficient 
of the beta function. Our default choice in Fig. [2a] is [3^ = IOO/Sq (/^o = 11/4 — nj/b) 
which is numerically close to the Fade estimate obtained in Ref. [3l|. For /34 = 0 and 
(^4 = 200/9o one observes a shift of the five-loop result by about -1-0.5 MeV and —0.5 MeV, 
respectively. 

It is interesting to look at the shift on at the central scale pdec = The two-, 

three- and four-loop curves lead to shifts of about —201 MeV, —21 MeV and —2 MeV, 
respectively. For = 100/3o the five-loop result leads to a shift of about —0.5 MeV. 

In a second application we consider the evaluation of {Mz) with (3 GeV) as input. 
The calculation proceeds in analogy to the bottom quark case discussed before, where for 
the on-shell bottom quark mass we use the value Mb = 4.7 GeV. Our results are shown 
in Fig. [2bl Again one observes a flattening of the curves after including higher order 
corrections. However, for p^®® 1 GeV, which corresponds to the left border of Fig. [2bl 

all curves show a strong variation which indicates the breakdown of perturbation theory 
for small scales. Around p^®®/Mf, ^ 0.3 both the four- and five-loop curves are basically 
flat. 

At the central scale p^®® = Mb one observes shifts in mf\Mz) of —55 MeV, —7 MeV and 
—1 MeV after including two-, three- and four-loop running accompanied by one-, two- 
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(a) 



(b) 

Figure 2: mf\Mt) as a function of (a) and m'c\Mz) as a function of (b). The 
numbers indicate the loop order used for the running. 


and three-loop decoupling. The shift at hve loops is below 1 MeV for /94 = 100/3o but also 
for /34 = 0 and [3^ = 200/3o- 
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3 Low-energy theorem: Higgs-fermion coupling 


The effective Lagrangian describing the coupling of a Higgs boson to gluons and light 
quarks can be written in the form 

£.([ = (Cio; + . (8) 

where the effective operators, which are constructed from light degrees of freedom [35], 
are given by 

0[ = , 

ni 

O'. = ( 9 ) 


The residual dependence on the mass rrih of the heavy quark h is contained in the co¬ 
efficient functions and C^- In Eq. ([H|) H denotes the Higgs held and v the vacuum 
expectation value. The superscript “0” reminds us that the corresponding quantities are 
bare. For the renormalization of (7°, C*^, O'^ and we refer to Ref. [II1|35]; for the pur¬ 
pose of this paper it is of no further relevance. In Ref. m a low-energy theorem has 
been derived which relates the computation of the renormalized coefhcient function C 2 to 
derivatives of (m w.r.t. the heavy mass mh- It is given by 


C 2 


^InCm 
d In rrih 


( 10 ) 


It should be stressed that Eq. fITOD is valid to all orders in ag- Note that Eq. ina contains 
the derivative w.r.t. Inrrih and furthermore the mh dependence of C 2 appears in the 
form ln{fi/mh)- Thus we can exploit renormalization group techniques to construct all 
logarithmic terms of the next, not computed perturbative order. In particular, on the 
basis of our four-loop calculation for (m we can compute C 2 to hve-loop accuracy using 
the recently computed hve-loop result for the quark mass anomalous dimension [8|. Note 
that the four-loop anomalous dimensions have been computed in Refs. [SI HI (7m) and 
Refs. OE] (/3), respectively. 

Inserting into Eq. (fTOj) we obtain the following result 
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where we have chosen /i = nih to obtain more compact expressions. Analytic result valid 
for general p are provided from [30] . 

In practice, one often encounters the situation where C 2 has to be inserted in a formula 
expressed in terms of If we furthermore transform the heavy quark mass to the 

on-shell scheme we obtain for /i = Mh 




0.2778 + 



3 

(1.355 + 0.2454 ni) 



4 

(-12.13 + 1.004 nz - 0.01003 



(-140.9 + 44.20 ni - 4.332 n^i + 0.07556 nf) . 


( 12 ) 


Let us briefly discuss the influence of C 2 on the Higgs boson decay to bottom quarks where 
the role of the heavy quark is taken over by the top quark. We consider the contributions 
proportional to {C 2 Y from Eq. ([8]) and use the result for the massless correlator from 
Ref. [M]- For convenience we identify the renormalization scale with the Higgs boson 
mass and set /i = Mh- Then the decay rate of the Standard Model Higgs boson to 
bottom quarks can be written in the form 

C 

V{H^hh) = -^^ml{MH)R{MH), (13) 

R{Mh) = 1 + 5.667 (^^^)+ (29.147 + 0.991) (41.758 + 13.105) (^^)^ 

+ (-825.7 + 50.7) + (rg + 224.8) (^^)^ (14) 

= 1 + 0.20400 + (0.03777 + 0.00128) + (0.00195 + 0.00061) 

+ (-0.00139 + 0.00009) + (O.OOOOOOOOrg + 0.00001) , 

with as = asiMn) ~ 0.1131. The first number in the round brackets in Eq. flT4)) cor¬ 
responds to the case C 2 = 1 [36] and the second one to the contribution from (C 2 — 1). 
At three-loop order the top quark induced part amounts to about 30%, at order aj only 
6%. Note that the massless correlator at order af, denoted by in Eq. ffT4|) . is currently 
unknown. The af term in Eq. fll4p origins from the hve-loop contribution in Eq. fll2p and 
products of lower-order contributions. 

Note that in this consideration the contribution of Ci (cf. Eq. (IH|)) has been neglected. 
The corresponding corrections of order can be found in Ref. EH. Corrections of order 
which are proportional to C'iC 2 require the evaluation of massless four-loop two-point 
functions and are currently unknown. Corrections of order to the Higgs boson decay 
rate involving (Ci)^ have been computed in Ref. [38] . 
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In Refs. [91110] the five-loop result for C\ is given in terms of and the MS quark 
mass. We complement this result by Ci parametrized in terms of the effective coupling 
constant and the on-shell mass: 


COS 


1 

12 TT 


1 + 



2.750 + 



(9.642 


0.6979n;) 



3 

(50.54 - 6.801ni 


0.2207n^) + 



- 625.2 + 149.8n, 


- 3.090nf - 0.07762nf + 6 1 , 


(15) 


where /i = Mf has been chosen. The analytic version in computer-readable form can again 
be found in [50] . 


4 Summary and conclusions 

In this paper we compute the four-loop corrections to the decoupling constant for light 
quark masses, Cm, which has to be applied every time heavy quark thresholds are crossed. 
It constitutes a fundamental constant of QCD and accompanies the hve-loop quark anoma¬ 
lous dimension [S] in the “running and decoupling” procedure. Our results complete the 
calculation of the four-loop decoupling constants which has been started in Refs. imn]. 
Note that the hve-loop corrections to the QCD beta function, which is needed to establish 
relations between Q;s(p) and mg(p) at low and high energy scales, is still missing. 

As a by-product of our calculation we obtain the effective coupling of a scalar Higgs boson 
and light quarks to hve-loop order. It is obtained from with the help of an all-order 
low-energy theorem. We briehy investigate the inhuence on T{H —)■ bb). 
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